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ABSTRACT 

The convective overshoot mixing plays an important role in stellar structure and evolution. However, the 
overshoot mixing is a long standing problem. The uncertainty of the overshoot mixing is one of the most 
uncertain factors in stellar physics. As it is well known, the convective and overshoot mixing is determined 
by the radial chemical component flux. In this paper, a local model of the radial chemical component flux is 
established based on the hydrodynamic equations and some model assumptions. The model is tested in stellar 
models. The main conclusions are as follows, (i) The local model shows that the convective and overshoot 
mixing could be regarded as a diffusion process, and the diffusion coefficient for different chemical element is 
the same. However, if the non-local terms, i.e., the turbulent convective transport of radial chemical component 
flux, are taken into account, the diffusion coefficient for each chemical element should be in general different, 
(ii) The diffusion coefficient of convective / overshoot mixing shows different behaviors in convection zone 
and in overshoot region because the characteristic length scale of the mixing is large in the convection zone and 
small in the overshoot region. The overshoot mixing should be regarded as a weak mixing process, (iii) The 
result of the diffusion coefficient of mixing is tested in stellar models. It is found that a single choice of our 
central mixing parameter leads to consistent results for a solar convective envelope model as well as for core 
convection models of stars with mass from 2M to 10M. 

Subject headings: convection — turbulence — stars: abundances — stars: interiors 



1. INTRODUCTION 

Many literature have shown that the convective overshoot 
mixing plays an important role in stellar structure and evo- 
lution. The classical treatment of the ov ershoot is based on 
the no n -local mixing le n gth theories (e.g..lShaviv & Sa lpeter 
(ll973l) : lMaederl (ll975l) : lBressan et al. I (119811) '). The convec- 
tive boundary in the framework of local mixing length theory 
(MLT) is the location where the acceleration of the fluid ele- 
ments is zero. In the framework of non-local mixing length 
theories, although the acceleration of fluid elements is zero 
(a = 0) at the convective boundary, the velocity is not zero, 
thus there should be an convective overshoot region outside 
the convective boundary. The non-local mixing length theo- 
ries trace the fluid elements overshooting from the convective 
unstable zone into the convective stable zone, find the loca- 
tion where the velocity of fluid elements is zero (v = 0), and 
regard the region laying between a = and v = as the 
overshoot region. The convective mixing in the overshoot re- 
gion, i.e., the overshoot mixing, is in usual assumed t o be very 
efficie nt and leads to complete mixing. However, iRenzini I 
dl987l) has argued that the theories of the mixing length type 
don't have enough spatial resolution to accurately describe the 
overshoot process, and the obtained overshoot distance in the 
mixing length theories is sensitive to the assumed turbulent 
heat flux in the overshoot region: V s» V ' a d leads to exten- 
sive overshoot distance and V ~ V/j leads to a small over- 
shoot distance. The problem can be resolved only with self- 
consistent nonlocal convection theory, in which the turbulent 
heat flux is not external imposed, but rather determined by the 



theory dGrossman et al. 1119931) . 

Another way to describe the overshoot mixing appeared 
in literature i s to regard the oversh o ot mixing as a d i ffusiv e 
process (e.g.. iFreytag et al. I (119961): iSalasnich et al.l (1999); 
lLai & Li I d201 lb ; Ehang & Li I d2012al) ). Thev set the diffu- 
sion coefficient of the overshoot mixing to be proportional to 
the local pressure scale height Hp and the turbulent veloc- 
ity. The turbulent velocity is set as an exponential decreasing 
functio n in the overshoot r egion according to numerical simu- 
~1 |199_6 ) and turbulent convection models 



lations (Frevtag et al. 
dXiong||1985L 



1989; Zhang & Li 2012b). Their representa- 
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tion of the adopted diffusion coefficient implies an assumption 
that the characteristic length of the overshoot mixing is on the 
order of Hp, which is same to the characteristic length of con- 
vective mixing in the unstable zone. However, this needs to 
be clarified. 

A more reasonable way to resolve the overshoot mixing 
problem is to use self-consistent nonlocal convection theory in 
multi-component fluid. The convective / overshoot mixing is 

determined by the turbulent chemical component flux u r 'X^ 
(u r is the radial turbulent velocity, Xk is the abundance of fcth 
chemical element, and u r 'Xk' is the cross-c orrelatio n func- 
tion of turbulent fluctuati ons of u r and Xk). IXiong I (11981b 
and ICanutol dl999tl20TTl) have developed the self-consistent 
nonlocal convection theories in two-component fluid for stel- 
lar interior based on some model assumptions. The domi- 
nating equation of the turbulent chemical component flux is 
established in those theories. Although this way is the most 
reasonable one among the three to deal with the overshoot 
mixing, it is also the most complicated one to be applied in 
the calculations in stellar evolution. More advances on this 
way are required in order to understand and apply it. 

In this paper, we present a model of the turbulent chemi- 
cal component flux u r 'Xk based on the hydrodynamic equa- 
tions and some model assumptions, and discuss the properties 
of the overshoot and the overshoot mixing. The derivation of 
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the model is presented in Section 2. The diffusion coefficient 
of the convective / overshoot mixing is presented in Section 
3. The property of the convective / overshoot mixing are dis- 
cussed in Section 4. The model is tested in stellar evolution in 
Section 5. Some discussions are presented in Section 6. The 
conclusions are summarized in Section 7. 

2. THE EQUATION OF THE CONVECTIVE FLUX OF CHEMICAL 
COMPOSITION 

Let us recall that the chemical composition thro ugh the fluid 
changes in two ways (lLandau & Lifshitz1ll987h : (i) caused 
by the macroscopic motion of the fluid, which leads to purely 
mechanical mixing, and (ii) caused by the diffusion, i.e., the 
molecular / atomic transfer of the components. In this pa- 
per, we focus on the first way, i.e., the purely mechanical 
mixing. In the absence of molecular / atomic diffusion, the 
composition of a fluid elem ent doesn't change when it moves 
(lLandau & Lifshitzll 19871) . namely: 



DXk 
Dt 







(1) 



where X k is the abundance of the fcth elements in the sys- 
tem (in this paper, 'k' means turbulent kinetic energy, but 'k' 
in subscript means the fcth chemical element), D / Dt is the 
Lagrange's derivation defined as follow: 
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where Ui describes the velocity vector of fluid. 

As it is usually used in the analysis of turbulence, a variable 
A is split into a mean and a fluctuating part as A = A + A' . 
We adopt Boussinesq approximation (i.e., taking into account 
the the density fluctuation only in buoyancy and ignoring 
the compressibility of fluids), which is generally accepted in 
studying stellar convection theories (IXiong et al. 1119971) . For 
this reason, this paper doesn't distinguish p and 7j, i.e., p = ~p 
in this paper. Split e quation (1) , note that ul — for the star in 
quasi-steady state (Xiong Tl989l) . and use Eq.(A7), one finds: 
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In the case of spherical symmetry, equation (1) becomes: 
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(pr 2 u r 'X k ') = 



(4) 



Equation (4) describes the convective mixing in stellar in- 
terior. It is clear that the convective mixing is determined by 

the radial chemical component flux u r 'X k . The goal of this 
section is to derive the representation of the radial chemical 
component flux. 
Let us start from the hydrodynamic equations: 
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where P is pressure, T is temperature, S is entropy, 
gi describes the gravitational acceleration vector, A = 
(AacT 3 ) I (3/c/o) is the thermal conduction coefficient, <Xjj 
describes the viscous stress tensor defined by Cjj = 
p(dui/ 'dxj+duj I 'dxi) + (£ — 2p,/ 3)5g(dui/dxi) in which a 
and £ are coefficients of viscosity (lLandau & Lifshitzl 1987). 
In the energy equation, T(DS/ Dt) can be replaced by: 



DS dS DT dS DP 

~ T{ df )pX ^t + T{ dP> T ' X ^t 
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where X — {X±, X2, X3, ...} is the chemical elements 
set including all independent elements (the sum of abun- 
dance of all chemical elements being unity indicates that 
there is one dependent element), Cp is specific heat, 6 = 
— (dlnp/dlnT) p t x is the expansion coefficient. Equation (1) 
is used in Eq.(8). The energy equation therefore can be re- 
written as: 



H Dt Dt 



du,. 



d_ 



dT 

dxj 



(9) 



In the Appendix, the evolutionary equation of second mo- 
ment is obtained, i.e, Eq.(AlO). Set A = Uj and B = X k (the 
abundance of an independent element) in Eq.(AlO), one finds 
the evolutionary equation of the convective chemical compo- 
nent flux: 
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where the mean Lagrange derivative is defined by D /Dt = 
(d/dt) + ul(d / dxi). According to Eq.(l) and (6), the above 
equation becomes: 
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The pressure fluctuation P' describes sonic process. For 
the subsonic convection, the sound wave is not important 
for transport pr ocesses (i.e., the anelastic approximation, see 
Spiegel ( 1971)) , thus we ignore the pressure fluctuation, fol- 
lowing [Xiong] (119811) . The density fluctuation in the buoy- 
ancy term is saved in above equation because the convection 
is driven by the buoyancy in stellar interior thus the buoyancy 
term p'gj is crucial. We use the formula as follows to model 
the density fluctuation in the buoyancy term: 



pf ,dlnp P> T' ^ 
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According to those, and note again Uj = 0, we can rewrite 
Eq.(ll)as: 
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Set A — T and B = X k in Eq.(AlO), one finds the evolu- 
tionary equation of the fluctuation T'X' k : 



DT'Xk' | l d(fmi'T'X k ') 
Dt p dxi 



(14) 
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According to Eq.(l) and (7), notes that the viscou s term 
is ignorable compa ring with the dissipation term (Kiong I 
198n!Canutoll997),and ignores the pressure fluctuation, one 
finds: 
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X k X k> in Eq.(13) can be modeled by analogy wifh lCamitol 
d20Tlh : 

WXV = ~C 1 T(^WX7+^WXV) (16) 
2 axi <JXi 

where r = k/e is time scale of the dissipation of turbulent ki- 
netic energy, k is turbulent kinetic energy, e is the dissipation 
rate of turbulent kinetic energy, C\ is a parameter. 
The dissipation terms are modeled as follows: 
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= C A (l + P e - 1 )r- 1 T'X k ' (17) 
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where P e = (pC p / \)(k 2 / e) is the Peclet number which de- 
scribes the ratio of the time scale of radiative thermal conduc- 
tion to the time scale of the dissipation of turbulent kinetic 
energy, Ca and Cp are parameters. 



In most cases of the stellar evolutionary stage, the stellar 
structure can be considered as in the quasi-steady state. In 
the quasi-steady state, the time derivation in Eq.(13) & (15) is 
zero. In additional, we assume that the turbulent fluctuations 
Uj'X' k and T'X k are in local equilibrium. Accordingly, the 

equations of the structure of Uj'X' k and T'X k in the quasi- 
steady local equilibrium (local model) are as follows: 
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+ (C B r- 1 + --^-C 1 cj ) T)u r 'X k ' 
z tip 

-^dCirYsVk'Ur'Xk/^ = 

fc' 

where = Vk'{dXy jd In P), g — —g r — GM r /r 2 

k' 

is the gravitational acceleration, V = dlnT/dlnP is the 
real temperature gradient in stellar interior, and V ' a d — 
(dlnT / dlnP) s is the adiabatic temperature gradient. 

3. THE DIFFUSION COEFFICIENT OF CONVECTION FNDUCED 
MIXING 

It is not difficult to obtain the the solution of the convective 
flux in Eq.(19) & (20) as follows: 



u r 'X k = -D k 



dXk 
dr 



(21) 



where 



D k = [C A (1 + Pe' L )T- L U r 'Ur' + i=U r 'T' (22) 

—CxCa^ + P e - 1 )gY / ur'Xk' 1 )] 

k' 

[C B C A (1 + P^T- 2 + \j-dCA{l + Pe-'W 



2 Hp 



Hp 



(V-V^)]" 1 



The value of Dk still can't be work out directly because the 

turbulent variable u r 'X k >' is present in the right hand side. 
However, the right hand side is independent of 'k' in the rep- 
resentation of D k - This indicates that Dk for all chemical 
element are the same in the local model, i.e., Dk = D. Ac- 
cordingly, one finds u r ' X k i' {dX k / dr) = u r ' X k {dX k * / dr) 
in Eq.(20) and then gets the representation of D: 

d = [ft(i + Pr 1 )r 1 iS' + =?r'] (23) 

{c B c A (i + Pe-'y- 2 

This result indicates that the convective / overshoot mixing 
can be regarded as a diffusion process in the local model. 
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The representation of Eq.(23) hints us to set C\ = 1/Ca- 
Note that the region with P e <C 1 is in the stellar surface 
convection zone and thus there is almost = 0, one finds the 
representation of the diffusion coefficient as follow: 



D = 



2ujC a (1 + Pe' 1 ) + 



Sg u r 'T' 



C A C B (1 + Pe- 1 ) + N 2 t 2 e 



(24) 



where to = u r 'u r ' / (2k) is the anisotropic degree, N 2 = 
~8g{y — V a d — 4>/5)/Hp describes Brunt- Vaisala frequency. 
The turbulent properties (i.e, the turbulent heat flux u r 'T', 
the turbulent kinetic energy k and the turbulent dissipation 
rate e) are required for calculating D. Those variables can 
be d etermined by the turbulent convection models ( TCMs) 
(e.g.JXiong)(fT^ : lCanutol (fT^:lDe" ng & Xiongl d20"06h: 
iLi & Yang I d2007l) : lLi I (120121) ). Here, we adopt Li & Yang's 
(2007) TCM. Set the model of uTT 7 and T T T ; to be local, 



i.e., Cti = and C e i = in the TCM, one obtains u r 
follow: 



Ur'T' 



2mC e {l + Pe- 1 ) 

CeCtil+Pe- 1 ) 2 +N 2 T 2 



-kr 



T' as 



(25) 



[-^(v-v ad )] 

where C e , Ct, C e \ and Ct\ are model parameters in Li & 
Yang's (2007) TCM (see Eq.(6), (7), (10) & (1 1) in their pa- 
per). 

Therefore, the diffusion coefficient of the convective / over- 
shoot mixing D is as follow: 



D = 2w(l + P e 



C A 



C c N 2 t 2 

C c C t (l + P c - 1 ) 2 + A' 2 T 2 



C B C A {l+Pe 



N 2 i 



(rk) (26) 



In stellar interior, the most important case for the convec- 
tive / overshoot mixing is in the region with P e ^> 1. Now, 
let us focus on the case of P e ^> 1. In the convection zone 
with P e > 1, Li & Yang's (2007) TCM shows N 2 t 2 
—C e C t /(2C e uj + 1) which can be worked out via simple al- 
gebra u sing the information in Appendix A in iZhang & Li I 
d2012bl) . In the overshoot region, TCMs show that the tem- 
peratu re gradient is close to the radiative t emperature gra- 
dient (IXiong I [l989t IXiong & Deng I l200lt iDeng & Xiongl 
l2006t rZhang & Li 2012b ), which is consistent with hel ioseis- 
mic investigation (IChristensen-Dalsgaard et al. 1 1201 ll) . thus 
N 2 t 2 ^> 1. Accordingly, a simple representation of D is 
as follows: 
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in the convection zone, and: 



D a 2lu(C a - C e ) 



N 2 



(27) 



(28) 



in the overshoot region. 



4. THE CHARACTERISTIC TIME SCALE AND LENGTH SCALE OF 
CONVECTIVE AND OVERSHOOT MIXING 

In the above section, it is found that the convective mixing 
can be deal with a diffusion process and the diffusion coef- 
ficient, i.e., Eq.(26), shows different behaviors in convection 
zone and overshoot region. In this section, we discuss the 



physical meaning of the diffusion coefficient and the mecha- 
nism of convective and overshoot mixing. 

In order to understand the diffusion mixing process, one 
should study the characteristic time scale tmix and the char- 
acteristic length of mixing Lyiix- The convection induced 
mixing is 'purely mechanic mixing', thus it is reasonable to 
regard the mixing as a consequence of the turbulent dissipa- 
tion. Therefore, the characteristic time scale of mixing should 
be comparable with the time scale of turbulent dissipation of 
kinetic energy, i.e., TMix ~ t. The characteristic length of the 
mixing should be the length of the radial range of the move- 
ment of fluid element, because the fluid element is finally dis- 
solved in this range. 

In the convection zone, the buoyancy impulses the radial 
movements of fluid elements, thus the radial length of the 
movement of an fluid element is the typical velocity multiplies 
the its lifetime which is comparable with the turbulent dissi- 
pation time scale of kinetic energy, i.e., Lmix ~ ykr. There- 
fore, the diffusion coefficient should be: D oc L 2 Mix /rMix ~ 
k 2 /e. 

In the overshoot region, fluid elements gain kinetic energy 
via turbulent diffusion of kinetic energy, and move (in ra- 
dial direction) around their equilibrium location, because the 
buoyancy prevents the radial movements of fluid elements. In 
this scene, the length of the radial range of the movement, also 
the characteristic length of the mixing Lmvx, can be estimated 
as Luix ~ y/k/N. Accordingly, the diffusion coefficient 
should be: D oc L 2 Mix /rMix ~ e/N 2 . 

It is necessary to discuss the relation between the turbu- 
lent heat transport and the matter mixing. In the high P e re- 
gion, the time scale of radiative thermal conduction is much 
longer than the time scale of the dissipation of turbulent ki- 
netic energy, thus there is almost no radiative heat exchange 
between fluid elements and the surrounding medium. There- 
fore, the heat exchange is a consequence of the matter mix- 
ing, or further, the turbulent dissipation. When the turbulent 
kinetic energy of a fluid element is completely dissipated, the 
fluid element is dissolved into the surrounding medium, thus 
it contributes its chemical composition and also its entropy to 
surrounding medium. The former causes the matter mixing, 
and the latter, i.e., entropy mixing, causes the heat exchange 
which is described by the turbulent heat transport. The tur- 
bulent mixing is not only on matter, but also on entropy. The 
turbulent heat transport can be regarded as a consequence of 
the matter mixing. When the mixing is effective, the entropy 
mixing leads to isentropic region, thus the temperature gradi- 
ent should be close to the adiabatic temperature gradient, and 
the chemical composition should be almost uniformed. On 
the contrary, when the mixing is ineffective, the temperature 
gradient should be close to the radiative temperature gradi- 
ent, and the chemical composition in the region, which is not 
chemically uniformed before taking into account the turbulent 
mixing, should be still not chemically uniformed. According 
to the discussions, the diffusion coefficient for heat transport 
should be similar to the case of matter mixing. Equation (25) 
shows that the diffusion coefficient Dt for the convective / 
overshoot heat transport is as follow: 



Dn 



2wC e (l + P e 



CeCtil+Pe- 1 ) 



N 2 i 



-(kr) 



(29) 



The representation of Dt also shows that Dt oc k 2 je in the 
convection zone and Dt oc e/N 2 in the overshoot region, as 
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similar as the case of matter mixing. 

Those phenomenological discussions explain the represen- 
tation of Eqs. (27-28). It is indicated in Eq.(26) that the 
convective mixing in the convection zone (|iV 2 r 2 | ~ 1) is 
much more efficient than the situation in the overshoot re- 
gion (N 2 t 2 ^> 1). This is because of the different char- 
acteristic length of mixing. In the overshoot region, the 
radial movements of fluid elements are around their equi- 
librium location, thus the characteristic length of mixing is 
much shorter than the situation in the convection zone. This 
could be the reason of why an excessively small parameter 
Cx = 10~ 10 is adopted in order to fit some observations 
when D = C x H P Vk (w C x k 2 /e, similar to Eq (27)) i s 
applied to the overshoot mixing (see Zhan g & Li I (l2012al) : 
iZhang I d20l2illbh \ 

5. APPLICATIONS IN STELLAR MODELS 

Usually, the turbulent kinetic energy k in the convection 
zone is on the magnitude order 10 6 or more, and the dissi- 
pation rate e = k 3 ' 2 /(aHp), where Hp is around the mag- 
nitude order 10 9 and a ~ 1. These lead to that the diffu- 
sion coefficient D is on the magnitude order 10 12 or more in 
the convection zone. The size of convection zone is about 
L ~ (1CP 1 ~ 10 1 )Hp. Therefore, the time scale of complete 
mixing in convection zone is about tMix ~ L 2 / D ~ (10 4 ~ 
10 s ) seconds or less, which is short enough in most cases in 
stellar evolution, leading to the complete mixed convection 
zone. The most interesting effect of the diffusion coefficient 
Eq.(26) is on the overshoot region. Equation (28) shows that 
the diffusion coefficient of overshoot mixing is proportional 
to [Ca — C e ), thus we define a new parameter: 

Cxi = C A -C e (30) 

In this section, we try to find the proper value of Cxi- 
The diffusive overshoot mixing based on the diffusion coeffi- 
cient Eq.(26) is tested in two cases: the solar model and the 
Hertzsprung-Russell (H-R) diagram of main sequence stars 
with convective core. The first corresponds to the convec- 
tive envelope downward overshoot mixing, and the second 
one co rresponds to the conv ectiv e core overshoo t mixing. Re- 
cently, Zhang & Li I (12012a!) and lZhangl (l2012 a.b) have stud- 
ied the overshoot mixing in stellar model of low-mass stars by 
using the diffusion coefficient described by D = CxHpy/k 
and found that the parameter Cx should be on the magni- 
tude order 10~ 10 . The excessively small parameter indicates 
that the representation of adopted diffusion coefficient, i.e., 
D = CxHpVk, seems to be not physically suitable. 

Li & Yang's (2 007) TCM is adopted in the s tellar modeling . 
iPaczvnskil dl969t) evolution code modified by Zh ang I d2012bl) 
is used to model stars. The modified code can solve the com- 
bination of stellar structure equations and the TCM in stellar 
evolution models. The diffusion coefficient of the overshoot 
mixing is describe by Eq.(26) i n the overshoot region. Param- 
eters of the TCM are same to Zhang (2012b), except for a. 
The value of a is 0.878 resulting from the solar calibration 
based on the TCM and the diffusion coefficient Eq.(26). The 
model parameter Cp is insensitive to the diffusion coefficient 
in the overshoot region. It is set to be Cr = "2Ct- 

The OPAL equation of state (iRogers et al. I 119961) . 
the OPAL opacity t ables for high temperatures 
(Igl esias & Rogers I 119961) . and the Alexander's opa city 
tables for low temperatures ([Alexander & Ferguso 1119941 ) are 
used. The composition mixture is assumed to be the same as 



the sola r mixture (|Grevesse & Sauval IH998I) . The elements 
settling dThoul et al. 111994b is used in the solar models. 

5.1. The convective envelope overshoot: the solar model 

For the solar model, there are two properties might be con- 
nected with the overshoot mixing below the base of the con- 
vective envelope, i.e., t he sound speed profile an d the surface 
Li abu n dance (see e.g..lSchlattl & Weiss I dl999t> ; iBrun et aTTI 
( 1999); Zhang & Li (2012a)). There are some different pro- 
posal to solve the Li pr oblem. However, the recent helioseis - 
mic investigation (i.e., iChristensen-Dalsgaard et al~l (1201 11) ) 
indicates that the convective overshoot mixing sho uld be 
taken into account when study the Li abundance dZhang I 
l2012al) . The solar models evolve from ZAMS to the present 
solar age 4.6Gyr. Three evolutionary series of the solar mod- 
els are calculated, include the standard solar model (MLT is 
used, no overshoot) and two solar models based on the TCM 
and overshoot mixing with Cxi = 10~ 3 and Cxi = 10~ 2 . 
The radius, luminosity and the boundary of the convection 
zone of the present solar models are calibrated to fit the ob- 
servations (R Q = 6.96 x 10 10 , L Q = 3.846 x 10 33 ) and 
helioseismic inversion (the base of the solar convection zon e 
R cz = 0.7135-Rm , IChristensen-Dalsgaard et al~l (fl99lh : 
iBasu & Antial d!997l) : iBasul d 1998b ) via the iterations of ini- 
tial hydrogen abundance, metallicity, and the TCM parameter 
a,. 

The overshoot mixing solar model with Cxi = 10~ 3 
shows 0.5 dex Li depletion in the main sequence stage, 
and the model with Cxi = 10~ 2 shows 1.4 dex. It has 
been f ound that the pre-main sequence Li depletion is about 
1 dex dD' Antona & Mazzitelli I [l994t iPiau & Turck-Chieze I 
2002; Zha ng & Li 1 12012a) . But the pre-main sequence Li 
depletion is se nsitive to the opacity, there is an uncertainty 
about 0. 5 dex dP' Antona &" Mazzitelli l 119941). The ob ser- 
vations dAnders & Grevesse II 19941: lAsplund et al. 112009b re- 
quire about 2 dex Li depletion. Accordingly, the main se- 
quence Li depletion should be 0.5 ~ 1.5 dex. This in- 
dicates that the parameter Cxi should be in the range of 
10~ 3 ~ 10~ 2 . The differences between the sound speed of 
the solar models and the helioseismic inversion are shown in 
Fig.l. It can be found in Fig.l that the sound speed has been 
improved in both two TCM solar models with overshoot mix- 
ing. The differences of the sound speed of model below the 
base of the convection zone are reduced to be less than 0.2% 
in both of the overshoot solar models. 

5.2. 77ie convective core overshoot: evolutionary tracks of 
main sequence stars in the H-R diagram 

For the stars with convective core, it is well known that the 
core overshoot mixing leads to significant effects on the stel- 
lar evolution. The classical treatment on the overshoot is to 
set an overshoot region outside the convective core and the 
the overshoot region is assumed to be complete mixed as well 
as the convective core. The length of the overshoot region 
is usually set as lov = &ovHp, where aov is a parame- 
ter. The proper range of the value of aov can be derived by 
compa ring the stellar models with corresponding observation 
data. iStothers I {1991) has found an upper limit of OAHp for 
the core oversh oot region (for the classically complete over- 
shoot mixing). Claret (2007) has investigated the length of 
the classical core overshoot region by using the observations 
of double-lined eclipsing binaries, and the results also indi- 
cated that, when the error bars are taken into account, the 
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FIG. 1. — The sound speed differences in the solar models. Sc 2 /c 2 = 
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FIG. 2. — Stellar evolutionary tracks in the HR diagram of in the main 
sequence stage. The dotted lines correspond to the standard stellar models, 
in which the convection heat flux is calculated by using the MLT and the 
overshoot mixing is absent. The dot-dashed lines correspond to the classical 
overshoot stellar models, in which the convection heat flux is calculated by 
using the MLT and the overshoot leads to complete mixing in OAHp outside 
the boundary of the convective core. The solid and dashed lines correspond 
to the stellar model with the TCM and diffusive overshoot mixing. The value 
of Cxi for solid and dashed lines is 10 — 3 and 10 -2 , respectively. 

core overshoot region based on the complete mixing seems 
to be less than OAHp. Claret's (2007) results and some of 
Stothers's (1991) are based on the properties of stellar evo- 
lutionary tracks in the H-R diagram. Therefore, their results 
require that the satisfied stellar evolutionary track of the star 
with a convective core should locate between the evolutionary 
track of OAHp completely overshoot mixing stellar model 
and the evolutionary track of standard stellar model without 
overshoot. 

Figure 2 shows the evolutionary tracks of 2M, 3M, AM, 
5M, 7M, and 10M stars. Stellar models in four cases 
are calculated: the standard stellar models (MLT is used, 
no overshoot), the classical overshoot stellar models (MLT 
is used, completely overshoot mixing in OAHp outside the 



Schwarzchild convective boundary), the TCM stellar models 
with the diffusive overshoot mixing with Cxi — 10 -3 and 
Cxi = 10~ 2 . The initial chemical composition is hydrogen 
abundance Xh = 0.7 and metallicity Z — 0.02 for all stel- 
lar models. The dotted lines correspond to the standard stel- 
lar models. The dot-dashed lines correspond to the classical 
overshoot stellar models. The solid and dashed lines corre- 
spond to the TCM stellar model with the diffusive overshoot 
mixing, where Cxi = 10~ 3 for solid and Cxi = 10~ 2 for 
dashed lines respectively. It can be found that only the the 
TCM stellar models with Cxi = 10~ 3 locate between the 
evolutionary tracks of OAHp completely overshoot mixing 
stellar models and the evolutionary tracks of standard stellar 
model. All of the TCM stellar model with Cxi = 10" 2 show 
too strong efficiency of overshoot mixing. Therefore, based 
on Li & Yang's (2007) TCM and adopted TCM parameters, 
the proper value of Cxi should be about 10~ 3 or less. 

6. DISCUSSIONS 
6.1. On the overshoot 



The ' ballistic' oversh o ot models (e.g. , IShaviv & Salpeter 
(119731) : iMaederl (119751) : iBressan et akl (119811) ) investigate 
how far a fluid element from the convective unstable zone 
can go into the stable zone, and define the location where the 
velocity of the fluid element is zero as the boundary of the 
overshoot region. In this description, the turbulent convective 
transport of turbulent kinetic energy is not taken into account. 
The overshoot flows stir the surrounding medium and trans- 
port turbulent momentum and kinetic energy, thus the medium 
originally locating in the stable zone should move. And, the 
moving medium go on stirring and transporting turbulent mo- 
mentum and kinetic energy. According to this mechanism, the 
real overshoot region should be larger than the prediction by 
'ballistic' overshoot models. This is consistent with turbulent 
convection models, which show very extensive overshoot re- 
gion that is much more extensive than the prediction by 'bal- 
listic' overshoot models. In this sense, one should regard the 
'overshoot region' as the overshoot of turbulent kinetic en- 
ergy rather than the overshoot of the fluid elements from the 
convective unstable zone. 

According to above discussions, the flows in the overshoot 
region can be classified as two kinds: (i) the flow overshoot- 
ing from the convective unstable zone ,and (ii) the flow caused 
by the turbulent convective transport of turbulent kinetic en- 
ergy. The 'ballistic' overshoot models focus only on the first 
type of flow. Fluid elements in the second type of flow are 
originally staying and locating in the stable zone. They gain 
turbulent kinetic energy via the turbulent convective transport. 
Because the buoyancy in the stable region prevents the radial 
movement of fluid elements, they should move (in the radial 
direction) around their equilibrium location. The radial range 
of their movements is much shorter than that in the convective 
unstable zone, thus the effect of the mixing (heat and matter 
transport) is much lower than that in the convective unstable 
zone. This indicates that, in the overshoot region, the tem- 
perature gradient is more closer to the radiative temperature 
gradient than to the adiabatic one, and the chemical composi- 
tion may not be uniform. The fluid element moves around its 
equilibrium location, thus the radial velocity is the maximum 
when the fluid element is in its equilibrium location where 
the turbulent temperature fluctuation is zero, and the radial 
velocity is zero when the turbulent temperature fluctuation is 
the maximum. This indicates that the correlation coefficient 
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between turbulent radial velocity and turbulent temperature 
fluctuation is much smaller than that in the convective unsta- 
ble zone. It is for the reason that, in the convective unstable 
zone, the fluid element is accelerated by the buoyancy and 
then turbulent temperature fluctuation and turbulent radial ve- 
locity increase in the same time, leading to significant correla- 
tivity between turbulent temperature fluctuation and turbulent 
radial velocity in the convective unstable zone. 

The low efficiency of the mixing (heat and matter trans- 
port) leads to V s» V# and weak matter mixing in the 
overshoot region. The former is consistent with turbulent 
convection models, which s how V »i Vp i n the overshoot 
region with P e > 1 (e.g., IXiong] (119851): IXiong & De"ng1 
d200ll) : lZhang'etal.1 (120121) : iZhangl (I2012bl) ). And, the 
latter is favored by the helioseismic results: the he lioseis- 
mic investigation (Christensen-Dalsgaard et al. 201 1) shows 
0.37 Hp o vershoo t below the base of the solar convection 
zone, and [Zhang (2012a) has found that complete mixing 
in the overshoot region with the length 0.37/Jp leads to 
excessive Li depletion in solar model. The low correla- 
tion coefficient between turbulent radial velocity and turbu- 
lent temperature fluctuation in the over shoot region is con - 
sistent w it h num erical simulations (e.g., Sin gh et aLl (1995); 
Meakin ( 2007 )) a n d turbulent convection models (e.g., 
Xiong I (11985lTl989l):lXiong & Chen I (fl992T) : lXiong & Deng! 



Along ( iy»y); Along &l Uien ( 1WZ); Along & 
(2001); Deng & Xiong I (120061) : iZhang & Li I d2012b)). 



6.2. On the parameter Cxi 

The results of Section 5 show that, based on Li & Yang's 
(2007) TCM and adopted TCM parameters, C X i = C A - C e 
seems to be about ~ 10~ 3 in order to satisfy some observa- 
tional restrictions. The value of Cxi being small indicates 
that Ca ~ C R . This is understandable. Let us recall the defi- 
nition of C P 



i.e. 



Eq.(7) in lLi& Yangl d2007): 



A dT' dT' 

pCp dxj dxj 



C e (l + P e - l )T- 1 T'T' 



(31) 



This is similar with the definition of Ca, i.e., Eq.(17). It 
should be noticed that AT (the difference of the tempera- 
ture between a fluid element and the surrounding medium) 
is strongly correlated with AXk (the difference of the abun- 
dance of fcth chemical element between the fluid element and 
the surrounding medium), since both of them are approxi- 
mately proportional to the radial distance between the location 
of the fluid element and its equilibrium location. This indi- 
cates that the turbulent temperature fluctuation is strongly cor- 
related with the turbulent abundance fluctuation, which means 
that T' is approximately proportional to X' k . Taking the strong 
correlation between T' and X' k into Eq.(31), comparing with 
Eq.(17), one finds Ca ~ C e . 

6.3. On the nonlocal model 

The diffusion coefficient in Section 3 is based on Eq.(19- 
20), which are the local limit of Eq.(13) & (15). An con- 
clusion based on the local model is that the diffusion coef- 
ficient for different chemical element is the same, namely, 
Dk = D. When the turbulent convective transport terms, 

i.e., {dpui'T'Xk'/dxi) and (dp Ul 'uj' X k ' / dx t ) in Eq.(13) 
& (15), are taken into account, this conclusion doesn't hold 
in general condition. This is because the turbulent convec- 
tive transport terms depend on high order radial derivations 
of chemical abundance. Therefore, in general, the efficient of 



convective and overshoot mixing is different for each chemi- 
cal element. However, there is a special case, i.e., when the 
abundance of ith chemical element Xi plus the abundance of 
jth chemical element Xj is a constant in the stellar interior. 
Xi and Xj satisfy dXi/dr = —dXj/dr because Xj + Xj 
is constant. In this case, Di = Dj always holds whether the 
turbulent convective transport terms are taken into account or 
not. This is because the equation: 



= Ur'{Xi + Xj)' = Ur'X/ + U r 'Xj' 

dx,, „ dx, ,„ dXi 



(32) 



B "—- D^ - -(D 



Or 



Or 



Or 



mathematically requires Di = Dj. An example is the calcu- 
lation of standard stellar model of main sequence star without 
the settling and atomic diffusion, in which the metallicity is 
uniform in the stellar interior and then the sum of the hydro- 
gen abundance and the helium abundance is constant. 

7. CONCLUSIONS 

In this paper, we studied the convective overshoot mixing 
in stellar interior. The convective and overshoot mixing is 

determined by the radial chemical component flux u r 'X k . 
A local model of the radial chemical component flux is ob- 
tained based on the hydrodynamic equations and some model 
assumptions. The model is tested in stellar models. The main 
conclusions are as follows: 

1 . The local model of the radial chemical component flux 
shows that the convective and overshoot mixing could be re- 
garded as a diffusion process, and the diffusion coefficient for 
different chemical element is the same. However, if the non- 
local terms, i.e., (dpUi'T 1 X k / dxi) and (dpui'uj' X k Jdxi), 
are taken into account, the diffusion coefficient for each chem- 
ical element should be in general different. 

2. The diffusion coefficient of convective / overshoot mix- 
ing shows different behaviors in convection zone and in over- 
shoot region. In the convection zone, the characteristic length 
scale of the mixing is large, thus the diffusion coefficient is 
high enough to ensure that the convection zone is completely 
mixed in most cases. However, in the overshoot region, the 
characteristic length scale of the mixing is much shorter than 
in the case of convection zone, thus the diffusion coefficient 
is much lower. The overshoot mixing should be regarded as a 
weak mixing process. 

3. The representation of the diffusion coefficient of mixing 
is tested in stellar models. It is found that, based on Li & 
Yang's (2007) TCM and adopted TCM parameters, the key 
parameter of overshoot mixing Cxi = Ca — C e should be 
about 10~ 3 in order to satisfy some observational restrictions. 
The small value of Cxi means Ca ~ C e , which is thought 
to be reasonable according to the properties of the convection 
motion in stellar interior. 
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APPENDIX 

SOME PROPERTIES OF THE FLUCTUATIONS 



Some basic properties: 



A = A + A' (Al) 



AB = A-B + A'B' (A2) 



A' = A-A = A-A = (A3) 



A'B' = A'B' + = A'B' + A'B = A'B (A4) 



(AB)' = (A + A')(B + B') - A- B - A'B' (A5) 
= A'B + B'A + {A'B 1 ) 1 

According to the continuity equation and Boussinesq approximation (ignoring the density fluctuation), one finds: 

d(pui)' 



As a consequence: 



The mean Lagrange's derivation: 



dx. 

, OA _ 1 d(pui'A) 
1 dxi p dxi 



= (A6) 



(A7) 



DA _ dA + _dA _ DA ^,dA 
Dt dt ' ' dxi Dt 1 dxi 
_ DA 1 d(pu/A) 



Dt p dxi 

The evolutionary equation of the second-moment: 



D(A'B') _ D(A'B') 
Dt ~ Dt 



(A9) 



AI DB' ni DA' DB DA 

A' h B'—— = A' h B' 

Dt Dt Dt Dt 



- A' DB A' ,dB B' DA B' ,dA 

Dt 1 dxi Dt 1 dxi 







+ A' Ul '— 

oxi 















-dB j^dA _ 1 d(pu t 'A'B') 



-A'm'- B 



dxi dxi p dxi 

In the derivation of the above equation, Eqs.(A4),(A7) & (A8) are used. The final result is as follow: 



D(A'B') l d(pu t 'A'B') 
Dt p oxi 
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The above equation describes the evolution of the second moment A'B' . The second term in the left hand side is the turbulent 
convective transport, the first two terms in the right hand side are the turbulent local productive terms, the last two terms are the 
contributions from the mean field. 
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